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Abstract 

The electric dipole response of the halo nucleus 11 Li is calculated in a hyper- 
spherical three-body formulation, and is studied as a function of the interaction 
employed for n— 9 Li to reflect the Pauli principle. Strength concentrations at lower 
energies are found but no narrow resonances. Only one possible scenario of 11 Li 
structure is in close correspondence with MSU and RIKEN experimental data. 
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Very recently, the experimental data for the electric dipole strength distribution dB(£l, E) / dE 
have become available from MSU JTJ for the two-neutron halo nucleus 11 Li at excitation 
energies E < 2 MeV. The Michigan fragmentation experiment for 11 Li + Pb was carried 
out at 28A MeV. While previous data on inclusive variables such as geometrical proper- 
ties and momentum distributions for individual break-up fragments have proven unable 
to discriminate between quite different wave functions for the bound state of n Li [Q, |3|, ffl, 
we will show that the dipole response function sheds some light on this issue. We remind 
the reader about the peculiar 3-body Borromean structure of n Li: since there is only 
one bound state (weakly bound at about 300 keV) and none of the binary subsystems 
are able to form bound states, the asymptotic conditions of the problem are pure 3-body 
asymptotics. 

Fig. 1 shows the deconvoluted experimental Michigan data which peaks at around 
0.6 MeV, with energy being measured from the break-up threshold. The same peak 
position was reported ||] for the RIKEN experiment with the same reaction at 42 A MeV. 
Fig. 1 also shows previous theoretical results, that of the simple point-dineutron cluster 
model of Bertulani and Bauer (BB)|| |7j, and the more realistic model of Esbensen and 
Bertsch (ESB)|8|]. We notice that the calculations, which nearly coincide, are shifted to 
low energies compared with the data. In the BB cluster model the energy dependence 
is simply E^J^E^I 2 j[E + E s ) 4 , where E s « 0.3 MeV is the separation energy of the two 
halo neutrons. This distribution peaks at E = |i? s , giving a linear dependence on the 
separation energy. The Green's function method of ESB employs a 5-type nn interaction. 

This letter investigates three limiting cases of plausible dynamics for the halo neutrons 
of 11 Li, corresponding to three tentative prescriptions for treating the Pauli Principle 
within a strict 3-body formulation. Expansion of hyperspherical harmonics (HH) is used 
||, a natural procedure for this Borromean system. 

We have recently JTIJ IT, [H| and with considerable success, carried out within the same 
HH framework bound state and continuum calculations for 6 He, also a Borromean nucleus, 
which shares many properties with n Li. Contrary to n Li, however, sufficient information 
is available on the binary ra-core interaction, and this makes physically reliable calculations 
possible. Both the + ground state and the known 2 + resonance at an excitation energy of 
1.8 MeV were reproduced by our calculations. However, although these calculations gave 
strength concentrations at low continuum energies, they showed no narrow resonances 
neither in the phase shifts nor in the calculated strength functions for electric dipole 
excitations of 6 He. A parallel study of n Li is hampered by lack of information on the 
n— 9 Li channel, although calculations can technically be carried out||. 
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The dipole response has the form 



£(£1;0+, - =| a(aL-(E!)\\7f ||0 g + s .) | 2 , 

where Tf M = I]j ejrjY" 1A f(fj) is the corresponding intrinsic 3-body cluster dipole operator 
with coordinates referred to the cm. of 11 Li, and hence free from spurious cm. motion. 
Note that the operator explores both the 1~ structure of the continuum as well as the 
structure of the + ground state. (We leave the assumed 9 Li(3/2 _ ) core out of all our 
calculations and discussion.) With previous experience from a similar calculation for 6 He 
|T2]| , we expect that replacing six-dimensional plane waves in the n — n and (nn)— 9 Li 
coordinates by correlated states from strict 3-body continuum calculations will shift the 
continuum strength to lower energies. Fig. 2 shows that this expectation is borne out 
in all cases. Notice that leaving out continuum correlations for the ESB case moves the 
peak into the position of the data, but that it now falls significantly below the data in 
strength. 

We follow [0, [13| and consider three choices, referred to as the H, WS-P and WS-F 



scenarios. For all cases considered, the nn interaction is that of Gogny-Pires-de Tourreil 
(GPT)|14| and includes repulsion at small distances, as well as spin-orbit and tensor 
forces. The Pauli principle is included by introducing repulsive 'Pauli cores' from the 
required combination(s) of central and spin-orbit forces. As discussed in previous papers 
H, U, all scenarios give a binding energy for n Li of ~ 0.3 MeV, i.e. within experimental 
error, and give [[| very similar single-particle densities. As in the 6 He case, investigations 
of the three-body scattering amplitudes fl5| showed no narrow dipole resonances states 



for any of the scenarios. These three scenarios give however very different ground-state 
wave functions, as exhibited in table 1. 

In the WS-F (Woods-Saxon-Forbidden) scenario, the spin-orbit force is assumed to 
play a leading role. We assume that there is a Opi/2 resonance in n + 9 Li scattering 
(|[§ 0, @), with the 0p 3/2 existing as a bound state of 9 Li (B.E.= 4.1 MeV). The 9 Li 
core is then taken to have a full level of neutrons, and the Pauli principle was taken 
into account by pure repulsive potential in the occupied 0s and 0p 3 / 2 states of neutron- 
core motion. The resulting 11 Li wave function (WF) will be predominantly (Opi/2) 2 , which 
corresponds (for harmonic oscillator orbitals) to a linear combination of 33% 1 S'o and 67% 
3 Pi states of neutron-neutron motion. This implies that the 11 Li ground state has a rather 
small probability of a di-neutron configuration, starting at 33% and increasing to 43% (see 
table 1) when correlations are included. 

In the WS-P (with P for pairing) scenario, we follow |1{J by taking the pairing rather 
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than the spin-orbit forces as dominant in n Li. A Cohen-Kurath type of calculation for 
9 Li gives |0| a ground state that is 93% of the symmetry [/] = [432], which can couple 
only with the spatially symmetric ([/] = [2]) neutron pair to form n Li with a closed shell 
structure. The valence neutrons are now almost entirely in (l nn — 0, S — 0) relative 
motion, with the S = 1 configuration blocked by the core neutrons. The spin-orbit force 
being damped in the interior so that it does not destroy the coherence of the state which is 
approximately y / 2/3(0p3/ 2 ) 2 + -y/l/3(0pi/ 2 ) 2 . We include this effect in our model by having 
no neutron-core spin-orbit forces, taking into account only the S = configurations in 
the Op shell. As in WS-F case we use a Pauli repulsive s-wave core-n interaction. 

The third (H) scenario has no explicit Pauli terms, but simply employs the shallow 



Gaussian potential of p0| |!|, which does not support any occupied orbits. Thus the mean 
fields for the halo and core neutrons differ substantially. The halo neutrons may now be 
largely in 0s motion relative to the core, and only partially blocked by Pauli orthogonality 
with the core neutrons because the core and valence radii are different. This calculation 
also reproduces the binding energy and r.m.s. matter radius of the n Li ground state and 
gives the largest value of hyper-radial moments <| p 4 |> compared to <| p 2 |> 2 (see Table 
1) as an expression of the softness of the n Li halo system. 

In all scenarios, the ground state of 11 Li has a closed neutron shell: (Os) 2 in the H case, 
and (0s) 2 (0p) 6 in the WS-P and WS-F cases. Thus, to form a 1~ excitation, a neutron 
will have to be excited into the next shell of opposite parity, implying an energy which is 
hu ~ 4 MeV in a standard estimate for dipole excitations. 

For halo nuclei, however, the neutrons are very near threshold, and simple estimates 
based on hu ~ 4 MeV may no longer be correct. Recently, for example, it was argued in a 
self-consistent density-functional method and also in a two-body cluster calculation 
[0, that the centre of gravity of strength functions should be concentrated closer to 
zero energy as the valence level approaches threshold. To answer this question, we have 
calculated continuum distributions in a 3-body model, which should include the effect of 
such shifts as a 3-body threshold feature. 

In all scenarios, we obtain a satisfactory convergence for the SI strength in the energy 
region below 10 MeV . However, for energies greater than about 4 or 5 MeV (the thresh- 
old for 9 Li excitation) more complicated mechanisms come into play, and our analysis 
assuming an inert core is no longer physically valid. We must also ensure convergence in 
the radial integrations. Since the effective range (a product of two wave functions and 
a radial operator) of the dipole operator is ~ 40 fm, to obtain the response without 
erroneous contributions from artificial low lying additional structures, we have extended 
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our calculations out to 150-200 fm. 

The results for continuum final states (using the same 3-body hamiltonian as for the 
g.s.) are shown in Fig. 1, compared with the MSU data. We show also the results 
(ESB) of the Green function method^], and the distributions (BB) of the 2-body cluster 
model|J. The WS-P and WS-F distributions, as expected from simple estimates, are 
very broad. The continuum RPA results |^2[ are similar to the WS-P curve. Only the H 
scenario resembles the MSU data. 

The curves peak at different energies. The position of the BB peak has already been 
discussed. In the WS-P model, qualitatively speaking, we have a purely repulsive s-wave 
n-core interaction, lifting up the levels in the next shell, and in WS-F there is also a 
repulsion, additionally pushing the levels to higher energies. The curves also approach 
E = in different ways. In the 2-body BB|7| point dineutron cluster model an E 3 ! 2 
behaviour is obtained for the low-energy SI response. Since the g.s. WF is concentrated 
in the asymptotic tail, it is possible to estimate this quantity in the three-body case, 
which gives ~ E 3 for the threshold behaviour. Both models give ~ E~ 5 / 2 when E is 
much larger than the binding energy. 

In the calculation of continuum WFs the ground state core-n potential and realistic 
GPT nn interaction were used. If we replace this nn interaction by a simple central 
Gaussian (which however reproduces n — n low energy phase shift data), all curves shown 
in figs. 1 & 2 are changed in the H case by less than the line width, and by only a 
few percent in the WS cases. Only the core-n and central part of the nn interaction are 
decisive for the n Li dynamics. 

In our calculations of the electric dipole response of 11 Li with 3-body wave functions 
for both the continuum and g.s., we summarise our results as follows: 

i) Although the strength is concentrated at lower energies, resembling resonant be- 
haviour, no narrow resonances are found in the three-body scattering amplitudes |T5| , 
implying a wide spreading of dipole states into the continuum; 

ii) Our calculations support qualitatively the conclusion that the enhancement of 
dipole strength at low energy is mainly due to the halo structure of the g.s. That is, 
the kinematical (threshold effect) enhancement is due to the large moments of the mass 
distribution, and to the proximity of the g.s. to the three-body core + n + n decay 
threshold. 

iii) 81 transitions which play the main role in electromagnetic dissociation are very 
sensitive to: 1) the structure of the ground state structure and to a lesser extent of 
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the continuum; 2) the method of treating the Pauli principle; and (not addressed here) 
possible mechanisms for the 11 Li dynamics during the reaction process. 

iv) Only the H-scenario is in a good qualitative agreement with the experimental SI 
deconvoluted response. It should be noted, however, that in all three cases the Pauli 
principle was taken into account in an approximate way. The question of how to handle 
the Pauli principle in halo nuclei is thus still open, and one we are still persuing. 

We will return to the nuclear monopole mode and the quadrupole excitation in a larger 
paper, where also sum rules are discussed in more detail. 

We wish to acknowledge support from NORDITA and the University of Bergen, where 
much of this work was carried out. U.K. support from SERC grants GR/H53648 and 
GR/H53556 is acknowledged. Helpful discussions with H. Schulz and experimental results 
communicated by him are also appreciated. 
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Figure 1: Electric dipole response for n Li from the MSU experiment [[J 
(squares), with the Greens function [|J (ESB, long dashed) and cluster model 
(BB, short dashed) calculations. The curves H (solid), WS-F (dotted), 
WS-P (dot-dashed) are the predictions from the current 3-body scenarios. 
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Figure 2: Electric dipole response with correlated (a) and uncorrelated (b) 
continuum wave functions, for n Li in the H, WS-F & WS-P scenarios, and 
the ESB calculation g. 
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Table 1: Partial norms, matter radii, binding energies for n Li g.s. in WS-P, WS-F and H 
scenarios, and mean energies for SI responses. K is the hypermoment quantum number, 
while l nn and l( nn ) c are orbital angular momenta in the n — n and (nn)— 9 Li degrees of 
freedom, and couple up to L. 
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